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Most of the electromagnetic output of blazars (BL Lac objects and Flat Spectrum Radio 
Quasars) comes out in the 7-ray band, making the Large Area Telescope [0.1-100 GeV] 
onboard the Fermi satellite and the Cherenkov telescopes crucial for gather crucial data 
and thus to understand their physics. These data are complemented by the observations 
of the Swift satellite in the X-ray and optical-UV bands, and by ground based optical 
and radio telescopes. This rich coverage of the spectrum allows a robust modelling, from 
which important trends start to emerge. In powerful sources we see the contribution of 
the accretion disk that, once modeled, give us the black hole mass and the accretion rate. 
Even when not directly visible, the disk luminosity can be derived through the broad 
emission lines. Therefore we start to know the jet power, the disk luminosity, and the 
black hole mass, 3 crucial ingredients if we want to draw a general scenario. At the start, 
jets are believed to be magnetically dominated. And yet, on the scale where they emit 
most of their luminosity, their power is already in the form of kinetic energy of particles. 
Relativistic jets are formed for a very broad range of the disk luminosity, from close to 
Eddington down to at least 10~ 4 Eddington. Their power correlates with the accretion 
rate, and can be even more powerful than the accretion disk luminosity. 
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1. Introduction 

The jets of Flat Spectrum Radio Quasars (FSRQs) and BL Lac objects, being rela- 
tivistic and pointing at us, are the most powerful persistent sources in the Universe. 
After decades of research we have still to explain the main mechanism launching, 
accelerating and collimating jets, but part of our partial lack of progress can be 
ascribed to our ignorance of basic physical parameters, such as the power of jets, 
their composition, the particle acceleration mechanisms and the efficiency in pro- 
ducing the radiation we see. In this respect, the Compton Gamma Ray Observatory 
( CGRO) with its high energy instrument, EGRET, and the ground based Cherenkov 
telescopes were crucial to unveil at least the entire apparent (i.e. isotropic equiv- 
alent) emitted luminosity, that often peaks in the 7-ray band. Now, with more 
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advanced Cherenkov telescopes and especially with the successor of CGRO, namely 
the Fermi satellite and its Large Area Telescope (LAT, 20 times more sensitive 
than EGRET), we entered in a new era (e.g. Ref. [TJ [3]). Besides the improved 
capabilities of these instruments, we also benefit from the flexibility of the Swift 
satellite, that makes possible to have routinely optical, near UV and X-ray data 
really simultaneously with the high energy ones. 

Furthermore, being visible up to large redshifts, blazars can be used as probes 
of the far Universe even if do not understand completely their physics. We can 
exploit the fact that the Universe is not transparent to high energy (~TeV) 7-rays, 
because of the cosmic IR background photons that can interact with 7-rays and 
absorb them, leaving some fingerprint of their presence. And, more recently, blazars 
open the search for the "cosmic magnetic field background" . 

Prospects are rather optimistic for blazar research, and in the following I will 
try to outline some of the results already obtained. 

2. Low power BL Lacs 

At the low end of the blazar sequence^^^ we find low power, TeV emitting, BL 
Lac objects. The TeV flux of PKS 2155-30^ and Mkn 501 6: showed, occasionally, 
ultrafast variability, with a doubling time of the order of 3-5 minutes. This, as 
discussed below, is a problem for our naive expectations, and challenges all models. 
Even more intriguing, and challenging, is the recent finding that not only lineless, 
low power BL Lacs vary fast, but even broad emitting line FSRQs 

TeV photons can be converted to pairs while traveling to Earth, by interacting 
with the cosmic optical-IR background photons. The resulting; absorption features 
have been used to test models of the optical-IR backgrounds, l2Qn23H14| Du t a morc 
recent use of this process concerns the estimate of the cosmic magnetic field. 

2.1. Ultrafast TeV variability 

Theoretically, one expects that the minimum variability timescale in jetted AGNs 
(and in Gamma Ray Bursts, too) is linked to the light crossing time of the 
Schwarzschild radius of the black hole, even if the emission region is very distant 
from the central engine. As an example, consider the "internal shock scenario" , 
popular in the Gamma Ray Burst field, but applicable to blazars too (for which, 
indeed, it was originally invented^) . Assume that two shells are initially separated 
by AR and have a similar width. Take AR of the order of the Schwarzschild radius 
Rs- Assume the front one moves with a bulk Lorentz factor T and the back one 
moves with 2r. They will collide at a distance i? C oii ~ T 2 AR ~ T 2 R$, and will 
interact for a distance equal to -R C oii- For observers on the jet axis, the Doppler 
contraction will shorten the observed time to i var = R co n/ (cT 2 ) ~ Rs/c. A similar 
timescale is expected for conical jets producing TeV photons at a distance ~ 10 3 i?s 
from the central engine, from a region with a size equal to the cross sectional radius 
of the jet. For a 10 9 M Q black hole, we expect t var ~ Rs/c ^3 hours. We instead see 
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hundreds of seconds, corresponding to a size R<ct val T ~3x 10 13 i V ar,2Pi cm. We 
then require very compact regions moving very fast. One possibility is to have small 
reconnecting regions, producing anisotropic currents of relativistic electrons in the 
comoving frame of the jet. We then have a "jet in a jet" scenarid^ (relativistic flows 
of particles in a frame that is also moving relativistically) . These electrons would 
inverse Compton scatter their own synchrotron radiation or any ambient photon 
field, cool very rapidly, and produce fast variations. An alternative model is to in- 
voke magneto-centrifugal acceleration of electrons along rigid and rotating B-field 
linea^ 21 . Electrons would reach their maximum speed close to the "normal" dissi- 
pation region, and could inverse Compton scatter very efficiently the synchrotron 
photons of that region. These electrons would move in a very collimated beam, and 
relativistic effect would be huge. A strong requirement of this model is that the 
accretion disk must be radiatively inefficient, since otherwise it would prevent the 
accelerating electrons to reach the required energies. In BL Lacs we do have ineffi- 
cient disks, and so one would be tempted to associate ultrafast TeV variability to 
BL Lac objects only. 

The recent detection^ of ultrafast variability (~10 minutes) from a FSRQ source 
by MAGIC (PKS 1222+216, z = 0.431) came therefore as a surprise. To let TeV 
photons survive, the emitting region should be away from the accretion disk (that 
in this source is efficient, we see strong broad lines), and yet we need a compact 
source. In this case electrons accelerated by the magneto-centrifugal force cannot 
reach energies high enough to produce TeV photons. On the other hand, also the 
"jet in a jet" reconnection model may have problems, since it requires a magnetically 
dominated jet, while the synchrotron luminosity of this source is not dominant (see 
Ref.|43]). 

2.2. Cosmic magnetic fields 

The TeV photons produced by blazars, absorbed by the IR background, produce 
pairs, that efficiently cool by Inverse Compton scattering off the CMB photons. 
Photons of energy e = 10 TeV would produce pairs of 5 TeV each, corresponding 
to a 7-factor of ~ 10 7 . Their reprocessed radiation will have e ~ 100 GeV. 

If the cosmic magnetic field is vanishingly small, the pairs will emit their radi- 
ation within the original emitting cone of the blazar. Therefore all the reprocessed 
radiation is contained in the same cone of the original TeV radiation. This implies 
that the flux is conserved: all the radiation absorbed in the TeV energy range must 
come out at lower frequencies. If, instead, there is a non-negligible i?-field, then 
the pairs would gyrate along the field lines while they cool, and the reprocessed 
radiation would be spread out in a cone larger than the original blazar's one. The 
received flux of the reprocessed radiation would then be lower. The left panels of 
Fig.Q]sketch these cases. Note that, in general, the high frequency reprocessed radi- 
ation, corresponding to high energy particles and/or the early phases of the particle 
cooling, is emitted in a short time, during which the particles do not gyrate much. 
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Fig. 1. Left: Cartoon of the reprocessing of the absorbed TeV radiation. TeV photons are 
emitted within a cone with semiaperture 8 C . They are converted into pairs at a typical distance of 
hundreds of Mpc. Pairs cool rapidly through inverse Compton scattering on the CMB. If B = 0, 
pairs do not gyrate, their emission is entirely contained in the original emission cone. In this case 
the observed reprocessed radiation has the same luminosity of the absorbed one. If B is non-zero, 
pairs with very high energy (top panel) cool before changing their direction (arrows) and thus their 
reprocessed emission is beamed within the same angle 8 C . Pairs of lower energies start to gyrate 
while cooling, and then emit within larger angles; thus their observed luminosity is smaller then the 
original one. The resulting overall spectrum can be approximated by three power laws: i) electrons 
emitting within the original cone (top), ii) within a cone larger than the original (middle) and 
iii) almost isotropically (bottom). Right: High energy SED of 1ES 0229+200. Lower high energy 
points are the observed H.E.S.S. spectrum, (Aharonian et al. 2007) higher points are de-absorbed. 
Lines are what predicted for different cosmic magnetic field values, as labelled, and for the initial 
cone angle 8 C = 0.1 rad and 0.05 rad. Black upper limits are from Fermi/LAT. From Tavecchio et 
al. 2010). 



As the particle energy decreases, the cooling timescale gets longer, the particle can 
gyrate more, and the reprocessed flux is spread within larger cones. When the gy- 
ration is complete, the reprocessed radiation becomes isotropic. We have therefore 
three regimes, corresponding to three different spectral slopes of the reprocessed 
flux, as shown in Fig. Q] (taken from Ref. |44|) . The right panel of the same figure il- 
lustrates an application to the TeV blazar 1ES 0229+200 (z — 0.14), a particularly 
interesting case because it has been observed up to ~ 10 TeV. The grey shaded 
region depicts the absorbed flux, and the different lines the expected level of the 
reprocessed flux for different cosmic i3-fields. The upper limits of the Fermi data 
limit the .B-field to be greater than ~ 10~ 15 G. As pionereed by Ref. l29l and 144} 
this technique is very promising, even if there are many complications that one has 
to account for with respect to the heuristic explanation given above (see e.g. Ref. 

mi). 
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Fig. 2. The luminosity of the BLR of blazars, in Eddington units, as a function of their 7-ray 
luminosity in the Fermi/hAT energy band, in Eddington units. A strong (linear) correlation is 
present. The grey horizontal stripe corresponds to iBLR/^Edd = 5 X 10 — 4 , and seems to divide 
(most of) BL Lacs from FSRQs. The upper limits on Z^blr/ ^Edd come from a sample of BL Lac 
objects present in the SDSS for which Plotkin ct al. (2011) could calculate the black hole mass 
through the M— a relation. They are all at 2 < 0.4 and show absorption lines. From Sbarrato et 
al. (2011). 



3. The FSRQ-BL Lac "divide" 

The "classic" and traditional way to classify blazars is through the equivalent width 
(EW) of their broad emission lines. FSRQs have EW larger than 5A (in the rest 
frame), while BL Lacs have lines with EW smaller than this value, or no lines at all 
(see e.g. Ref . HT|) . This criterion is simple and of immediate practical use, and it can 
be though of a sort of measure of the thermal to non-thermal ratio. On the other 
hand, we know that the non-thermal jet continuum is highly variable, and that 
most of the emission is at high energies, not in the optical, and therefore the EW of 
optical lines do not measure the ratio of the bolometric thermal and non-thermal 
luminosities. 

There are blazars classified as BL Lacs that have lines of small EW, but whose 
luminosities is even larger than the one found in FSRQs. Examples are PKS 0208- 
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51 JI1 A q 0235+164 34 , PKS 0426-380^1, PKS 0537-44ll3Il. Furthermore, the high 
amplitude optical variability makes some FSRQs appear as BL Lacs when their non- 
thermal continuum is in the high state, and conversely, BL Lacs in a particular faint 
state can show broad emission lines that, albeit weak, can have EW greater than 5 
A (as BL Lac itselfS3). 

There is the need of a more physically based classification scheme, and we^l 
have proposed a distinction based on the luminosity of the broad emission lines 
measured in Eddington units: Lblr/ -^Edd- Normalizing to the Eddington luminosity 
ensures the appropriate comparison among objects of different black hole masses. 
We proposed that when £BLR./£Edd<^5 x 10~ 4 the objects are FSRQs, and are 
BL Lacs below this value. The original sample of Ref. [24] was enlarged by Ref. [37] 
confirming the earlier result. Fig. [2] shows the result of the latter study. Note the 
large number of upper limits both in £blr/ -^Edd and in L 7 /LEdd- These refer to 
BL Lac objects selected by Ref. [32] to be rather nearby (z < 0.4), with signs of 
stellar absorption by stars of the host galaxy in their SDSS spectra. They allowed 
Ref. [35] to measure the black hole mass through the M-a relation. Note that they 
all lie in the "right" quadrant of Fig. [5] i.e. Lblr/ -^Edd < 5 x 10~ 4 . Fig. [5]shows a 
good correlation between Lblr and the 7-ray luminosity in the Fermi energy range 
(the correlation remains essentially the same also when using absolute values, not 
normalized to Eddington), and a continuity between BL Lacs and FSRQs. 



3.1. ADAF? 

This poses a question: do the subdivision between BL Lacs and FSRQs reflects a 
fundamental difference of only a different overall power (of both the jet and the 
accretion disk)? We have proposed earlier (Ref. [22] for Fermi blazars; Ref. [T9l for 
FRI and FRII radio-galaxies) that one key and fundamental difference could be 
the accretion regime: FSRQs (with strong emission lines) should have geometrically 
thin, optically thick accretion disks accreting above Ld/L^d — 10~ 2 , while BL 
Lacs could have disk accreting below this value, and have Advection Dominated 
Accretion Flows (ADAF; e.g. Ref. [28]). The found division at L BLR /L E dd ~ 5 ~ 10 -4 
is roughly consistent with Ld/^Edd ~ 10~ 2 , especially if the covering factor of the 
BLR is ^1/20. This idea fits with the overall spectral properties of blazars: strong 
and efficient disks ionize the broad line clouds, and the broad lines can be an efficient 
source of external seed photons to be scattered by the jet to form the dominating 
high energy hump. Weak and inefficient disks imply weak lines, fewer seed photons, 
less cooling for the emitting electrons, and a weaker high energy hump, produced 
mainly by the synchrotron self-Compton process, but peaking at higher frequencies, 
because the emitting electrons, suffering less cooling, can attain higher energies. 
Furthermore, as discussed more below, the power of the jet likely correlates with 
the mass accretion rate: weaker disk have weaker jets. 

On the other hand, the SED of ADAF is remarkably different from the "mul- 
ticolor" black-body of a standard^ disk, with a very reduced fraction of ionizing 
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photons^. And yet we do observe broad emission lines even in rather extreme BL 
Lacs, such as Mkn 501, Mkn 421 and PKS 2005-489 (see Ref. HI and references 
therein), thought to have heavy black holes and very weak disk. We do not have 
enough examples to draw any strong conclusion, but note that the ADAF hypoth- 
esis, explaining the BL Lac/FSRQ divide, is not mandatory. The transition to the 
ADAF regime could occur at much smaller values of Ld/^Edd, as suggested by e.g. 
Ref. 1411 One simple alternative is to assume that the location i?diss m the jet where 
most of the radiation is produced is roughly the same in all blazars, once measured 
in Schwarzschild radii (i.e. i?diss ~ 10 3 -Rs). This is what comes out from fitting their 

IOOI 

SEDs . The relation between the radius of the BLR and the ionizing luminosity 

1/2 i 

-Rblr oc i ion (i.e. Ref. [26j [T0|) means that a weak disk (even if not ADAF) implies 
a small -Rblr- If -Rblr becomes smaller than -Rdiss, the external radiation, as seen 
in the jet comoving frame, becomes negligible. In this case the division between 
BL Lacs and FSRQs is still dependent on Ld/^Edd but, in addition, there is also a 
dependence on the black hole m ass, making the expected divide more blurred (but 
still consistent with existing datsP^). 




4. Jet and accretion power 

Since the jet radiation is beamed, it is not trivial to estimate the jet power from 
what we see. In a pioneering studjBS, the average jet power was estimated from 
the extended radio emission, calculating the minimum (equipartition) energy of the 
radio lobes (or extended structures) , and then dividing by an estimate of the source 
lifetime. Very interestingly, the jet power was found to linearly correlate with the 
luminosity of the narrow emission lines (thought to be produced by photo-ionization 
by the accretion disk), suggesting an approximate equality of the jet power and the 
disk luminosity. At large scales, the jet power can be measured in optical (by HST) 
and/or X-rays (by Chandra). Another promising way to measure the jet power is to 
use the X-ray bubbles/voids seen in the intra-cluster medium hosting radio- AGNs, 
allowing to measure the PdV work, the age, and then the jet power 21 " 1 ^. 

However most results on the power of blazars' jet consider the zone of the jet 
where most of the radiation is produced, assuming that this zone is a single one (a 
blob), characterized by the same magnetic field, particle density and bulk motion. 
Through modelling, we can derive how many leptons and magnetic field are required 
to account for the observed emission, what is the size, and the bulk Lorentz factor. 
The jet power is in several forms, but all can be calculated through: 



where R is the size of the emitting blob, U[ is a comoving energy density. The 
different forms of power are the magnetic one (U' = Ub), the kinetic power of the 
emitting electrons (U' = U e = n e (7)m e c 2 ), the kinetic power of the cold protons 
(U' = Up = n p m p c 2 ), and the power in the produced radiation \U' — f7 ra d ~ 
L' hol / (AttR 2 c)]. To calculate the latter we need the observed bolometric luminosity 




Pi = TrR 2 T 2 pcU[ 



(1) 
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and an estimate of the bulk Lorentz factor: 

r I p2 t obs j obs 

P - nR 2 T 2 Bc bQl - bQl ~ (2) 

where the last equality assumes that the Doppler factor S = T. This is a lower limit 
to the total jet power Pj et , since the jet would stop if all its (kinetic+Poynting) 
power is used to produce the radiation we see (even if in extreme TeV BL Lacs this 
can be exactly what happens: this would explain why F is very large in the radiative 
region, but small on the larger VLBI scale as indicated by mildly superluminal, or 
even sub-luminal apparent speeds; see Ref. l20 | [T8|) . 




10 43 10 44 10 45 10 46 10 47 10 43 10 44 10 45 10 46 10 47 

L d [erg s-'] L d [erg s" 1 ] 



Fig. 3. The power spent by the jet to produce the radiation we see, P T (left) and the total 
jet power Pj et (right) as a function of the accretion disk luminosity Lj, calculated through the 
observed broad emission lines (assuming = IOLblr)- Pt can be considered as a robust lower 
limit to Pjeti which is here calculated assuming one cold proton per emitting electron. Grey solid 
lines indicate equality. 



Fig. [3] shows P r (left panel) and Pj C t (right) as a function of the accretion disk 
luminosity for a sample of Fermi blazars. It contains both FSRQs and BL Lacs. 
For FSRQs the accretion disk component can be often seen directly in the optical- 
UV range; for BL Lacs it is deduced mainly from the luminosity of the (weak) broad 
emission lines. Fig.[3]shows that: i) the jet power is large; ii) it nicely correlates with 
id! if one takes one proton per emitting electron (see Ref. HH 02] for a discussion of 
this issue), then P ]ct > id, on average, suggesting that 

F jct ~ Mc 2 (3) 

for all values of M . Sometimes, it can be even larger (it surely is during exceptional 
flares, as in 3C 454.3^2H11). The fact that P ]ct correlates with Ld, but that it can be 
larger than that, is an apparent paradox: the correlation would point to accretion 
as the prime mover of jets, but then how can Pj e t be larger than Ld? This paradox 
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can be solved by noting that, to extract the black hole rotational energy, we need 
a strong magnetic held, that has to be produced/sustained in the inner parts of 
the accretion disk. We then need B% ~ 8irpoC 2 in the innermost regions of the 
disk to make the Blandford-ZnajeH^ process efficient. This links the jet and the 
accretion powers. This is also borne out by recent numerical simulations (e.g. Ref. 
|46|) . showing that, occasionally, Pj e t can be even larger than Mc 2 . 

Table 1. Summary of the average physical properties of Fermi bright blazars. Values from 
Ghiscllini et al. (2010). Pmi ss /R s is the distance of the dissipation region from the central 
black hole of mass M (in Schwarzschild radii). T is the bulk Lorentz factor. B is the magnetic 
field in the emitting region, as measured in the comoving frame. is the luminosity of the 
accretion disk. P' is the power injected in the form of relativistic electrons, measured in the 
comoving frame. Pj et is the total power of the jet, assuming one proton per emitting electron. 
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5. Conclusions 

Tab. Q] summarizes the average physical properties of BL Lacs and FSRQs, derived 
from fitting their overall SED (including the 7-rays). The main difference is their 
overall power (both of their jets and their disks), inducing a difference in the value 
of the magnetic field in the jet emitting region. Since the black hole masses are 
similar, the difference of the accretion power means a difference in the Ld/^Edd 
ratio, and a well defined divide. The division between BL Lacs and FSRQs occurs 
at a given luminosity of the broad emission lines measured in Eddington units: 
Lblr/ ^Edd = 5 x 10~ 4 , which corresponds to Pd/^Edd = 10~ 2 if the covering 
factor of the BLR is ~ 1/20. The changing regime of accretion (from radiatively 
inefficient to efficient) can be the explanation of the different "look" of BL Lacs and 
FSRQs. Alternatively, the disk could remain radiatively efficient even for small 
values of Pd/pEdd, and in this case the BL Lac/FSRQ "divide" could be due to the 
emitting region being outside (BL Lac) or inside (FSRQs) the Broad Line Region, 
whose size scales as L 1 ^ 2 . 

A useful lower limit to the jet power is P r , the power spent to produce the radi- 
ation we see. It is a robust estimate because it depends on the observed bolometric 
luminosity and the bulk Lorentz factor, so it is nearly model-independent. Since the 
7-ray luminosity is often the dominant contribution to the bolometric one, instead 
of L b oi one can use L^, and the following relation: 

P r = 1.3 x 10 7 P° 8 ergs" 1 (4) 

where P 7 is the K-corrected 0.1-10 GeV luminosity. This relation is shown in Fig. 
|4]together with the observed points. It can be used to quickly estimate a lower limit 
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Fig. 4. The power P r that the jet spends to produce the radiation we see as a function of the 
7~ray luminosity in the Fermi/LAT energy band (K-corrected). The bisector of the least square 
best fit yields P r oc Z/JJ 8 (and L 7 oc P r 1,27 ). The non linearity is mainly due to the fact that for low 
power BL Lacs the 7-ray luminosity (as measured in the Fermi/LAT energy band underestimate 
the bolometric one. 



to the jet power, bearing in mind that i) the 7-ray luminosity can vary by even 2 
orders of magnitude and ii) that the real jet power Pj e t can much larger (~ 2 orders 
of magnitude) than P r , if charge neutrality is due to cold protons. 

Blazars (and then their misaligned counterparts) have jets for a large range 
of accretion disk luminosities, in Eddington units. They range from ~1 (close to 
Eddington) to ~ 10~ 4 . The SED of blazars is often dominated by the high energy 
hump, with the synchrotron luminosity being a minor part in FSRQs. This limits 
the value of the magnetic field. Indeed, it turns out that the Poynting flux in the 
emitting region is less than P r . This means that in FSRQs the jet power is mainly 
in the bulk motion of cold protons. 

On the other hand, the only option to accelerate jet in blazars is through the 
Poynting flux: contrary to Gamma Ray Bursts, the jets of blazars are initially 
optically thin. The scattering optical depth tt is less than unity. Even if e pairs 
were copiously produced, they would quickly cool and annihilate, making tt<1 
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at the very start. This means that any internal radiation cannot be transformed 
into bulk motion (this can instead occur in Gamma Ray Bursts). Internal energy 
in the form of very hot (i.e. 7>15) matter would again very quickly cool in the 
dense radiation field produced by the nearby accretion disk. Again, we cannot use 
this internal energy to accelerate blazars' jets. The only remaining option is then 
a strong initial Poynting flux, that must be dominant at the start, and yet it must 
become a minor component in the emitting region, ~ 10 3 Schwarzschild radii away. 
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